A New Measurement of the Radiative K e3 Branching Ratio and Photon Spectrum 
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present a new measurement of the branching ratio of the decay Kl 
with respect to Kl — > it e^v e (K e z), and the first study of the photon en- 
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ergy spectrum in this decay. We find BR(K c37 ,E; > 3OMeV,0 c % > 20°)/BR(K o3 ) = 
0.908 ± 0.008(stat.)^QQ 13 ,(syst.)%. Our measurement of the spectrum is consistent with inner 
bremsstrahlung as the only source of photons in K e 3~t . 



I. INTRODUCTION 

Precise measurements of the properties of radiative de- 
cay modes of kaons test theories of kaon structure. Addi- 
tionally, understanding the phenomenology of these de- 
cays enhances the ability to perform other precision mea- 
surements and rare decay searches which rely on robust 
background predictions of which radiative decays may 
be a component. There are two distinct components in 
most radiative decays: direct emission (DE), where the 
photon is irreducibly part of the decay interaction; and 
inner bremsstrahlung (IB), where the photon is emitted 
from an external charged leg. The IB component is well 
understood and dominant in most radiative decays. It 
is the size and structure of the DE component which is 
important for understanding kaon structure. 



Fearing, Fishbach and Smith (FFS)§, § |] and 
DoncelM performed the first theoretical studies of radia- 



tive K e3 decays (K L 
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e37 ) 



FFS present 



the matrix element with the IB component and a phe- 
nomenological model of the DE component: 
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where Tib is the IB portion of the matrix element, M is 
the kaon mass, e is photon polarization, I is the electron- 
neutrino current vector, k is the photon momentum, K is 
the kaon momentum and Q is the pion momentum. They 
present a calculation of the relative K^ branching ra- 
tio for photon energies above 30 MeV in the kaon center 
of momentum frame (CM), using current algebra to es- 



2 



Photon Veto Detectors 



20 cm I 



Muon Muon 
Analysis Magnet Filler Veto 

Csl 



K L I 

beams 



Vacuum Decay Region 



Vacuum Window 



Drill 

Chambers 



Trigger 
Houoscopes 

Hadron Veto 
with Lead Wall 



140 160 180 

Distance from Target (m) 



FIG. 1: The KTeV Detector, as configured for collecting the 
data used in this analysis. 



timate the size of the four parameters in Equation (Q). 
FFS estimate a DE component roughly 1% the size of the 
IB component, but do not include it in their prediction 
of the K e 3 7 branching ratio. Doncel performs a similar 
calculation of only the IB component, emphasizing in ad- 
dition the need to cut on the angle between the electron 
and photon in the CM in order to avoid a background 
from bremsstrahlung photons generated in the detector 
material H. 

More recently, HolsteinQ and Bijnens et a/.|| have 
performed Chiral Perturbation Theory (%PT) calcula- 
tions. Holstein finds a smaller DE component than FFS, 
by a factor of 5-10, but like FFS does not include it in a 
numerical prediction of the branching ratio. His predic- 
tion for the branching ratio is the same as FFS. Bijnens 
et al. do not explicitly separate their calculation into IB 
and DE pieces. 

The NA31 collaboration]^] performed the most re- 
cently published measurement of the K e ^ branching ra- 
tio. With approximately 1400 events, they measured 
a radiative fraction BR(K o37 ,E 7 > 30 MeV,t9* 7 > 
20°)/BR(K o3 ) = 0.934 ± 0.036±g;g||%. With over 15000 
K e 3y events, KTeV approaches the sensitivity required 



to see the DE component in this mode. 



II. THE KTEV DETECTOR 

This measurement was performed using the KTeV de- 
tector at Fermilab, shown in Figure § KTeV was de- 
signed to measure 3?(e'/e), the direct component of CP 
violation in neutral kaon decays, which requires excellent 
detection capabilities of both charged and neutral par- 
ticles. Hence, KTeV is an ideal apparatus for studying 
decays like which have both charged and neutral 

particles in the final state. The aspects of the apparatus 
pertinent to this analysis are described below. 

The Ki, beams used in this experiment were produced 
by an 800 GeV beam of protons striking a beryllium oxide 
target. Collimators and sweeping magnets downstream 



of the target produced two, nearly parallel neutral beams. 
The decay volume for accepted decays began 110 me- 
ters from the target to allow the K$ component to decay 
away, and continued to the first drift chamber of the spec- 
trometer at 160 meters. The acceptance for decays up- 
stream of 122 meters was restricted by the "Mask Anti" 
(MA) anti-coincidence counter which had holes to allow 
the beams to pass. The kaon beams and decay products 
traveled through vacuum from the target to the first drift 
chamber. 

The KTeV spectrometer measured the charged decays 
products. It consisted of four rectangular drift cham- 
bers, each with two horizontal and two vertical planes of 
sense wires, and a large dipole magnet which imparted a 
transverse momentum of 0.412 GeV/c. The drift cham- 
bers measured horizontal and vertical track position with 
a resolution of 110 /im and momentum with a resolution 
of 0.4% at a typical momentum of 36 GeV/c. 

Energy measurements and particle identification were 
performed by a pure cesium iodide electromagnetic 
calorimeter (Csl). It consisted of 3100 blocks in a square 
array 1.9 m on a side and 0.5 m deep. Two 15 cm 
square beam holes allowed the passage of the neutral 
beams through the calorimeter. The calorimeter was cal- 
ibrated using momentum-analyzed electrons, with aver- 
age energy resolution for electrons of 0.75% (where the 
momentum resolution has not been subtracted out). The 
calorimeter was also used to associate tracks between ver- 
tical and horizontal views. 

Following the Csl was a five meter long iron muon fil- 
ter. Muons needed to have a momentum greater than 7 
GeV/c in order to traverse the filter and deposit energy 
in the veto scintillation plane. 

A series of "photon veto" counters surrounded the fidu- 
cial volume of the detector to detect particles missing 
the Csl. These veto counters suppressed the background 
from Kl — > 7r + 7T~7r° (K n a) and other decays with extra 
photons. 

The data for this analysis were collected in a dedi- 
cated low intensity (10% of nominal e'/e intensity) run 
during a 24 hour period in 1997. The low intensity of 
these runs was crucial to having a manageable accidental 
background component to our measurement. Only min- 
imal trigger conditions were applied: trigger hodoscope 
signals indicating the presence of at least two charged 
tracks, in anti-coincidence with photon veto counters. 



III. ANALYSIS 

The analysis of this decay proceeded in two stages. 
First, an inclusive K e 3 sample was isolated, then a sub- 
sample containing one photon was identified. 
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FIG. 2: Data/MC comparisons of the total charged momenta 
(P en ). The x 2 / D OF = 95.5/72. 

A. Kez Sample 

The presence of an electron and a pion originating from 
a common vertex define the K e3 normalization sample. 
Particle identification was performed by comparing the 
momentum measured by the spectrometer with the en- 
ergy deposited in the calorimeter: the electron was re- 
quired to have 0.95 < E/P < 1.1; the pion, E/P < 0.7. 
Events with extra tracks were rejected. The momentum 
of each track was required to be above 7 GeV in order to 
reject muons that stop in the muon filter before reach- 
ing the muon veto counter. Fiducial cuts were applied 
to tracks and clusters to restrict particles to well under- 
stood parts of the detector. Activity in the photon veto 
counters was not allowed (deposited energy < 0.3 GeV), 
however, there was no cut on extra energy deposits in the 
Csl and photons could go down the beam hole without 
being detected. A plot of the charged particle momen- 
tum of the normalization K e % sample is shown in Figure ^ 
along with a Monte Carlo simulation of the same. 

1 The Quadratic Ambiguity 

We need to know the photon momentum in the CM for 
both the branching ratio and spectrum measurements. 
However, because the neutrino in K e3 decays is unob- 
served, events cannot be reconstructed unambiguously. 
Kinematic constraints allow one to determine the mag- 
nitude of the neutrino's CM momentum along the kaon 
flight direction but not the sign (forward or back). This 



leads to a two-fold ambiguity in determining the kaon's 
lab momentum, and hence the photon's CM momentum. 

The acceptances for the two kaon lab momentum solu- 
tions are different, favoring one solution over the other. 
We picked the more likely solution for any pair of mo- 
menta based on MC calculations of the acceptance. This 
procedure resulted in the correct solution 65% of the 
time. 

The quadratic ambiguity is smallest in cases where the 
neutrino emerges perpendicular to the kaon flight direc- 
tion. Because of the finite resolution of the spectrom- 
eter, cob 2 (9* Kv ) can have non-physical, negative values. 
However, events in which the kaons scatter in the colli- 
mators can also have negative values of cos 2 (9* Ku ). To 
retain most of the unambiguous events while rejecting 
large, unphysical values from events in which the kaon 
scattered, we required cos 2 (9* Kll ) > —0.2. 

2 Kinematic Cuts and Backgrounds 

To increase the purity of the K e % sample, we required 
P| < 0.055 GeV 2 and M en < 0.5 GeV, the kinematically 
allowed values. 

iVr3 events were a background at the tenth of a per- 
cent level in the K e3 sample because of the small prob- 
ability (0.3%) for a pion to be misidentified as an elec- 
tron. K V 3 events were more significant (a few percent) 
for the A" e 3 7 subsample because the decay also has pho- 
tons present. The mass of the unseen 7r° in decays, 
restricts the kinematics of the charged particles. By elim- 
inating events in which the unobserved 7r° would have a 
physically allowed momentum when analyzed as we 
removed 99% of this background, while removing only 
7% of the K e3l signal. 

Likewise, well over 99% of Kt, — > 7r + 7r~ events with a 
misidentified pion were removed by cuts on the two pion 
invariant mass (494-502 MeV) and transverse momen- 
tum (< 50 MeV 2 ). 

A few misreconstructed decays remained in the 
sample and were subtracted using a MC sample of K^s 
events. This sample was normalized to the data in the 
case where the kinematic requirement above has been 
dropped and the signal had been estimated with a sample 
of MC K e 3 events. A small K e 4 background was also 
subtracted using a MC sample normalized to the inferred 
data kaon flux. These two subtractions together change 
the normalization sample by less than 0.1%. 

B. Kez-, Subsample 

The K e 3y subsample was selected from events in the 
inclusive K e3 sample. The subsample consisted of events 
with exactly one photon candidate cluster. 

The photon candidate cluster was required to be more 
than 8 cm from the electron cluster in the Csl and more 
than 40 cm from the pion cluster. The photon-electron 
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FIG. 3: Data/MC comparisons of the photon energy in the 
laboratory frame at two different scales. The photon energy 
was required to be above 3 GeV to be included in the 
sample. This cut is indicated by the vertical line in each plot. 
The x 2 /DOF = G.2/U for the first plot and 24.4/27 for the 
second. 



distance cut allows for reliable separation of the electron 
and photon clusters. The photon-pion distance cut re- 
moved events in which clusters due to pion shower fluc- 
tuations, which are not well modeled in MC, might have 
mimicked a photon shower. The energy of a photon can- 
didate cluster was required to exceed 3 GeV, while the 
transverse profile of the cluster was required to be con- 
sistent with the hypothesis of a single, electromagnetic 
shower. A plot of the photon energy in the laboratory 
frame is shown in Figure |3|. The photon cluster was 
also required to be more than 2 cm away from the elec- 
tron's position as projected from the upstream track seg- 
ment to the Csl. This cut removed photons from physical 
bremsstrahlung as the electron passed through the detec- 
tor. 

The -ftT e 3 7 events were required to have exactly one pho- 
ton candidate cluster. In addition, the calorimeter was 
required to be free of other electromagnetic clusters with 
energy greater than I GeV and which were separated 
from the electron cluster by more than 4 cm and from 
the pion cluster by more than 30 cm. The MC does not 
simulate the multiplicity of these clusters perfectly. A 
small correction, CV, of 0.25% was applied to account 
for this. The correction was determined by finding the 
ratio of the number of events with two or more photon 
candidate clusters to the number of events with one or 
more photon candidate clusters. 

For the branching ratio measurement, additional kine- 
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matic cuts on the energy, E* > 30 MeV, and the angle 
between the electron and photon in the CM, #* 7 > 20°, 
were used in order to obtain a result with a kinematic 
range commensurate with the previous experimental re- 
sult and to emphasize the interesting kinematic region 
for DE studies. 

In our study of the photon CM spectrum (see Section 
VI), the cut on the distance between the photon candi- 
date and the electron projection was lowered to I cm, to 
allow us to use looser kinematic cuts E* > 25 MeV and 
> 5°. 

The background subtractions mentioned above for the 
inclusive K e ^ sample were more significant in the 
subsample, comprising 0.7% of the total. 



IV. MONTE CARLO SIMULATION 

The detector acceptance used in calculating the 
branching ratio and in determining the acceptance- 
corrected photon spectrum was calculated by MC sim- 
ulation. As mentioned above, the photon momentum in 
K e 3 7 cannot be constrained due to the unseen neutrino. 
This ambiguity makes precise MC simulation essential 
to the success of this analysis. To accurately simulate 
accidental activity in the detector, a special trigger was 
used to record events at random times with a frequency 
proportional to the overall beam intensity. These random 
events were overlaid on MC generated events for compar- 
isons with the data and for acceptance calculations. We 
were aided by the low beam intensity conditions under 
which the data were recorded. 

FFS|| present a full listing of all terms for the squared 
matrix element, including both IB and DE components. 
We used their listing explicitly in our Monte Carlo (MC) 
simulations. For the branching ratio measurement, all 
the DE coefficients were set to zero. The MC sam- 
ple was combined with a MC sample of non-radiative K e ^ 
decays at a level corresponding to NA31's branching ra- 
tio measurement Q . The cutoff for generating a physical 
photon to trace through the detector was I keV in the 
CM. Photons below this energy as well as loop effects 
were taken into account in the K e ^ generation. The K e ^ 
form factor, used for both radiative and non-radiative 
events was A + = 0.0274(14) ||, based on a preliminary 
KTeV measurement. Our result depends only weakly on 
the form factor, and the form factor used is in agreement 
with the 2000 PDG value §. 

Two other features of the MC simulation are worth 
mentioning. First, electromagnetic and hadronic showers 
in the Csl were handled by separate libraries of show- 
ers generated with GEANT. Second, the tracking of 
particles through the detector and beamline accounted 
for a number of effects including multiple scattering, 
bremsstrahlung and synchrotron radiation, delta parti- 
cle emission, photon conversion and kaon regeneration. 
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Normalization, K e 3 


RawEvents 


5760888 


Ktt3 background 


267 


K e 4 background 


197 


Backgroundsubtractedevents 


5760424 


K e 3 Acceptance 


0.1067 


Signal, K e3l {E* > 30 MeV, 


e* ei > 20°) 


RawEvents 


15575 


Ktt3 background 


50 


K e 4 background 


62 


Backgroundsubtractedevents 


15463 


K e 3-, Acceptance 


0.03161(14) 



TABLE I: Sample sizes used in the branching ratio calcula- 
tion. 

V. BRANCHING RATIO CALCULATION 

The measured relative branching ratio is simply de- 
fined by the ratio of the events we measure in each sam- 
ple corrected by the acceptance in each case. It is impor- 
tant to realize that our dependence on MC is mitigated 
by the ratio of acceptances appearing in the calculation. 
The sample sizes and acceptances are listed in Table |. 
Thus we find, 

T(K e3l , E* > 30 MeV, e * 7 > 20°) 

= N(K e37 ) Acc(K e3 ) 

N(K c3 ) Acc(K e37 ) V 
= 0.908 ± 0.008(stat)±g;^(syst)% ( 2 ) 

As mentioned above, the correction, Cy = 1.0025 ± 
0.0009, is necessary to take into account the different 
rates at which K e3l candidate events are vetoed by the 
cluster multiplicity cut in data and MC. Vetoed events 
in the MC are not included in the MC calculation of the 
acceptance and thus require an explicit correction. 

Our result is compared with the other published mea- 
surements and predictions in Figure ^|. The uncertainties 
for the theoretical predictions are based solely on the ac- 
curacy of the stated results. FFS give a result only for 
6** 7 > 0. We have used our MC (based explicitly on the 
FFS H matrix element) to calculate the fraction of the 
events with #* 7 > 20° and give the result for this case. 

A. Systematic Uncertainties 

Estimates of the sizes of various significant systematic 
uncertainties are listed in Table |J. The largest system- 
atic uncertainty comes from the variation of the calcu- 
lated branching ratio with the E 7 cut, as the cut was 
varied by ±1 GeV about 3 GeV. Photons below 2 GeV 
are not well reconstructed as we have seen by analyzing 
large samples of decays. The relative size of this 
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FIG. 4: A comparison of our branching ratio measurement 
with other recent measurements and predictions. 

uncertainty reflects the steeply falling photon spectrum 
at low energies combined with the difficulty of modeling 
the calorimeter at the lowest energies. 

The error associated with picking the right momentum 
solution was based on the degree of consistency between 
the subset of the data where momentum solution is un- 
ambiguous within resolution, and the data set as a whole. 
The error on the acceptance was based on the statistics of 
MC generated. To this was added the variation in accep- 
tance due to the unknown DE parameters studied below. 
The upper limit on the background due to clusters origi- 
nating from pion shower fluctuations was determined by 
studying decays in data. 

The contribution of accidental activity in the Csl to 
the K e3l subsample was small, 0.5%, as determined by 
generating MC with and without accidental overlays and 
comparing the number of events with exactly one photon 
candidate. The number of events with more than one 
photon candidate agreed between data and MC with ac- 
cidental overlays at the 20% level, whereas there were dif- 
ferences of orders of magnitude without accidental over- 
lays. 

Other errors were estimated from the differences be- 
tween data and MC shapes of various distributions com- 
bined with different expected shapes between normaliza- 
tion and signal samples, and the uncertainty in the mea- 
sured K e 3 form factor. 
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Down(%) 


u P (%) 


E 1 cut 


-0.005 


+0.008 


K momentumsolution 


-0.007 


+0.007 


MCStatistics 


-0.006 


+0.006 


PionshowerBG 


-0.003 


+0.000 


K e z formfactor 


-0.003 


+0.003 


P* shape 


-0.003 


+0.003 


shape 


-0.002 


+0.002 


VetoCorrection 


-0.001 


+0.001 


PhotonConversions 


-0.001 


+0.001 


cos 2 (8* Kv ) shape 


-0.001 


+0.001 


AccidentalOverlay 


—0.001 


+0.001 


K W 3 BGsubtraction 


-0.001 


+0.001 


K e 4 BGsubtraction 


-0.001 


+0.001 


TotalSystematicUncert ainty 


-0.012 


+0.013 


StatisticalUncert ainty 


-0.008 


+0.008 


Tot alU ncert ainty 


-0.014 


+0.015 



TABLE II: A list of significant uncertainties in the branching 
ratio measurement. 

VI. SPECTRUM MEASUREMENT 

With the statistics available in this experiment, it 
becomes feasible to compare the FFS model for direct 
emission (Eq. (|l|)) with our observed photon CM mo- 
mentum spectrum. We do not have the sensitivity to 
consider all four parameters, so we use the soft-kaon 
approximation |l(J to set A = B = and restricted our 
study to C and D whose terms remain non-zero in this 
approximation. The soft-kaon approximation assumes 
that terms in the decay matrix element proportional to 
the kaon rest mass are negligible. 

To determine the values of C and D favored by our 
data, we used MC to generate spectra at a set of points 
in CD space. We compared these spectra, bin-by-bin, 
to the acceptance-corrected spectrum seen in the data 
to calculate a % 2 value at each point. The only free 
parameter in this comparison was the normalization of 
the two spectra. To increase the statistics for these fits 
the photon-electron angle cut was relaxed to Q* ei > 5°. 
This cut is still far from allowing physical bremsstrahlung 
events into our sample. The fit was done for photon en- 
ergies between 25 and 200 MeV in seven bins. Below 25 
MeV we were limited by the sensitivity of our calorimeter 
to soft photons. Above 200 MeV, the error due to making 
the wrong kaon momentum choice became dominant. 

The constant x 2 contours corresponding to a = 1 
(A 2 - Xmin = 1) and cr = 2 (x 2 - xLnJ= 4) as a func- 
tion of C and D are shown in Figure A polynomial 
interpolation was done to estimate the values x 2 between 
sampled points. We find C = — 5 ± 10 and D = 5+2i> 
with a strong correlation between the two parameters. 

We are most sensitive to the linear combination of C 
and D which approximates a linear perturbation of the 
IB-only spectrum. The precision of our result is more 
apparent when one chooses axes aligned with this combi- 
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FIG. 5: The o = 1 ( X 2 - Xmin = 1) and a = 2 ( X 2 - Xmin = 4) 
contours of the fit/interpolated \ surface as a function of C 
and D. The dot shows the position of the minimum, while the 
cross shows the IB-only point. The region within the oblique 
rectangle is rotated and redisplayed in the inset. The new 
axes C and D' are obtained from C and D after a rotation 
of 6 = 25.8°. 



nation. Choosing axes C and D' rotated by 8 — 25.8°, 
we find C = -2.5±Hand D' = 6± 22 . This is illustrated 
in the inset of Figure a. 

The observed acceptance-corrected E* spectrum is 
shown in comparison with the generated spectrum at the 
best-fit lattice point in Figure ||. The x 2 value for this 
point is 4.3 with four degrees of freedom. At the IB-only 
point (0,0), the x 2 is 7.4. 

Knowledge of the acceptance correction is the domi- 
nant systematic error contributing to our measurement 
of C and D' . Physical backgrounds are very low and in- 
significant in the spectrum measurement. We estimated 
the systematic uncertainty associated with C (the large 
statistical uncertainty in D' makes the systematic uncer- 
tainty insignificant) by multiplying the measured accep- 
tance profile by a linear factor and renormalizing to keep 
the overall acceptance constant. The allowed size of the 
linear factor was determined by statistical consistency 
with the measured acceptance. By varying the accep- 
tance in this manner we find a systematic uncertainty in 
C of 1.5 and in D' of 1.0. 



VII. CONCLUSION 

We have measured the relative branching ratio 
for radiative K e3 decays and fit the photon spec- 
trum to place limits on the size of DE components 
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FIG. 6: Data/MC comparison of the acceptance-corrected E* 
spectrum at the best CD point. The x 2 /DOF = 4.3/4, where 
the 0-25 MeV bin has not been included in the comparison. 

in the K e ^ matrix clement. Our measurement, 
BR(K o37 , E; > 30 MeV, 0* 7 > 20°)/BR(K c3 ) = 0.908 ± 

0.008(stat)^o]oi2( s y s t)%i is nearly five times more precise 
than the previous measurement, and in agreement with 
it. It is significantly lower than all published theoretical 
predictions. 

This is the first attempt to measure DE terms by 
studying the photon spectrum. The spectrum measure- 
ment is consistent with IB as the only source of pho- 
tons in if e37 decays. In the soft-kaon approximation, 



we find C = —5 ± 10 and D = 5_ 2 ]\ or > m a frame 
rotated by 25.8°, C = -2.5t^(stat) ± 1.5(syst) and 
D' = 6±24 ± l.O(syst). 
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